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Abstract The possible associations between the microquasars LS 5039 and 
LS I +61 303 and the EGRET sources 3EG J1824-1514 and 3EG J0241+6103 
suggest that microquasars could also be sources of high-energy 7-rays. In this work, 
we present a detailed numerical inverse Compton (IC) model, based on a micro- 
quasar scenario, that reproduces the high-energy 7-ray spectra and variability 
observed by EGRET for the mentioned sources. Our model considers a population 
of relativistic electrons entrained in a cylindrical inhomogeneous jet that interact 
through IC scattering with both the radiation and the magnetic fields. 
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Microquasars are a selected class of X-ray binaries that produce relativistic radio jets (Mirabel 
& Rodriguez 119991 Fender |2004). The origin of the jets is related to the matter accreted by the 
compact object, a neutron star or a black hole, from the companion star. These systems behave 
as scaled-down versions of quasars and active galactic nuclei. The population of microquasars 
is still a very reduced one, with about sixteen known objects up to now (Ribo 2002). Due to 
the presence in these systems of both a relativistic jet with a population of highly relativistic 
leptons and strong radiation fields, microquasars are promising 7-ray emitter candidates through 
IC interaction (for synchrotron self Compton (SSC) models, sec, i.e. Atoyan & Aharonian 
1999; for external Compton (EC) models, see, i.e. Kaufman Bernado et al. I2002J) . although 
these sources could emit at such energy ranges by other means (i.e. see the hadronic model of 
Romero et al. I2TTU3) . We are interested on determining whether a microquasar jet emitting by 
IC scattering can reproduce the EGRET data of the two sources associated to microquasars: 
3EG J1824-1514/LS 5039 and 3EG J0241+6103/LS I +61 303. Our numerical model considers 
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a population of relativistic electrons entrained in a cylindrical inhomogeneous jet interacting 
with both the radiation and the magnetic fields, accounting for the synchrotron, the EC and 
the SSC electron energy losses. 

2 LS 5039 AND LS I +61 303 

LS 5039 is a High Mass X-ray binary system (HMXB) located at 2.9 kpc from Earth (Ribo 
ct al. 2002J) . The stellar companion is a bright (V ~ 11) star of ON6.5 V((f)) spectral type, the 
compact object seems to be a neutron star, the orbital period of the system is P — 4.4267 ± 
0.0005 days, the eccentricity, e = 0.48 ± 0.06, and the semi-major axis of the orbit a = 2.6 x 
10 12 cm (McSwain et al. l2()()4|l . The microquasar nature of LS 5039 was clearly established when 
non-thermal radiation produced in a mildly relativistic jet was detected, being also proposed 
as the counterpart of the unidentified 7-ray source 3EG J1824— 1514 (Paredes et al. [2000). 

LS I +61 303 is another HMXB whose optical counterpart is a bright (V~10.8) star of 
B0 V spectral type (Paredes & Figueras ll986f) . presenting and equatorial disk (Be star) and 
harbouring likely a neutron star as a compact object (Hutchings & Crampton I1981f) . This 
source is at a distance of about 2 kpc. The more accurate value for the orbital period is 
P = 26.4960 ± 0.0028 days (Gregory EEEJ), the eccentricity is e - 0.7 (Casares et al. l2lTOj) . 
and the orbital semi-major axis is about a = 5 x 10 12 cm. This source presents variability at 
different timescales, from radio to 7-rays (see, i.e., Taylor & Gregory F1982I Paredes & Figueras 
HM)1 Goldoni & Mereehetti ITMl Tavani et al. HMSI Massi l2"UtI3|) . The microquasar nature of 
LS I +61 303 was established when a mildly relativistic jet structure was detected through VLBI 
observations of this source (Massi et al. 120011 12004(1 . althought this object had already been 
proposed to be the counterpart of the 7-ray source 2CG 135+01/3EG J0241+6103 (Gregory & 
Taylor E75I Kniffen et al. lTM7|> . 

3 MODELLING THE 7-RAY EMISSION FROM MICROQUASARS 

In this model, we assume that the leptons of the jet dominate the radiative processes related to 
the 7-ray production. The relativistic population of electrons, already accelerated and flowing 
away into the jet, is exposed to external photons as well as to the synchrotron photons emitted 
by the electrons, since we take into account the magnetic field in our model. The 7-ray emitting 
region, the 7-jet, is assumed to be closer to the compact object than the observed radio jets. 
This 7-jet is supposed to be short enough to be considered cylindrical. The magnetic field (-B 7 ) 
has been taken to be constant, as an average along the jet. The energy losses of the relativistic 
leptonic plasma at the 7-jet are mainly due to synchrotron emission, SSC scattering, and EC 
scattering. Due to the importance of the losses, the electron energy distribution density along 
the 7-jet model varies significantly, thus the 7-jet is studied by splitting it into cylindrical 
transverse cuts or slices. 

Regarding IC interaction, we have used the cross-section of Blumenthal & Gould l|197()(l . 
cr(x, eo, 7c), which takes into account the Thomson and the Klein-Nishina regimes of interaction; 
eo is the seed photon energy, 7 C is the scattering electron Lorentz factor, and x is actually a 
function which depends on both of the former quantities and on the scattered photon energy 
(e). The electron distribution is assumed to be initially a power law (N(j c ) oc E~~ p , where E 
is the electron energy), which evolves under the conditions imposed by the magnetic and the 
radiation fields. Thus, the electron distribution function of a certain slice (N(^ c ,z)) depends 
on both the distance to the compact object (z) and 7 C . The components of the total seed 
photon radiation field {U{eo,z)) are any present external radiation field (Z7ext(eo; Z Y) an( i the 
synchrotron radiation field produced by the jet's relativistic electrons, all of them in the jet's 
reference frame (for the external photon fields, see Dermer & Schlickeiser I2()02jl . It must be 
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noted that, since LS 5039 and LS I +61 303 are not strong X-ray emitters, we will assume that 
the most important source of seed photons is the companion star, neglecting at the present stage 
the disk and corona contributions (see Bosch- Ramon et al. 12004(1 . The formulae that represent 
£/(eo, z) and its components (a black-body stellar photon field and a synchrotron photon field) 
are well described in the work of Bosch- Ramon & Paredes 120041 

The free parameters of the model are B 1 and the maximum electron Lorentz factor at 
the slice closest to the compact object (7^o ax )- The leptonic kinetic luminosity or leptonic jet 
power (ike) is set free also, and it is scaled with the observed luminosity, in order to reproduce 
the observations. The jet radius, although it is not well constrained, has been fixed to several 
electron Larmor radii (~ 10 7 cm), and the electron power-law index of the injected distribution 
has been deduced from radio observations (for LS 5039, see Marti et al. 1998 for LS I +61 303, 
see Ray et al. 11997(1 . The luminosity per energy unit in the reference frame of the jet (L e ) 
is presented in Eq. ^ The photon flux per energy unit or spectral photon distribution in the 
reference frame of the observer (I' e ,) is shown in Eq. [21 The magnitudes with (') are in the 
observer reference frame. 

L t = e > Kiice(z) / / N{j Cl z) d-yde (1) 



r < = I^v L *' (2) 

where 8 is the Doppler factor of the jet, D is the distance from the microquasar to the observer, 
and V s n cc (z) is the volume of the slice at a distance z from the compact object. For further 
details of the model, see Bosch- Ramon & Paredes I|2()()4ll2?)04b| . 



4 RESULTS AND CONCLUSIONS 

We have applied our model to both LS 5039 and LS I +61 303 in order to reproduce the 
spectra observed by EGRET. In Tabled the known parameter values as well as the adopted 
values for the free parameters are presented. The computed spectral photon distributions are 
shown in the Figs. ^ an d HI for LS 5039 and LS I +61 303, respectively. The results show 
that these microquasars can produce high-energy 7-rays through IC interaction between a 
relativistic leptonic population, entrained in an inhomogeneous cylindrical jet, and the external 
and synchrotron photon fields. 

Our model can reproduce EGRET data with reasonable parameter values, and variability 
is naturally expected from stellar radiation and accretion changes due to orbital eccentricity 
(precession can introduce also variability, see Massi et al. 12004(1 . Since, at the present stage, 
no clear varying accretion influence has been detected in LS 5039 at lower energies, we have 
not introduced it to compute the presented spectral photon distributions (see Fig.Q}, but only 
the eccentric orbit effects on the stellar radiation density. However, since orbital eccentricity 
is clearly affecting radio (Taylor & Gregory [1982|) and X-ray emission (Goldoni & Mereghetti 
1995) in the case of LS I +61 303, we have taken into account it to calculate its emission at 7- 
rays. The magnetic field can control the loss of energy of the electrons by SSC losses. Therefore, 
with two magnetic field strenghts: 1 and 10 G, the emission turns from being dominated by 
EC effect to be dominated the SSC effect. It is shown in Fig. ^ where emission is not affected 
by the eccentricity of the orbit when the magnetic field is high, since changes in the stellar 
photon density do not affect the spectrum in the SSC dominant regime (accretion changes are 
not included in the calculations). It is worth noting also that in the case of LS I +61 303, where 
accretion changes are taken into account, they are far more important than the stellar photon 
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Table 1 Parameter values for LS 5039 and LS I +61 303 



Parameter 


Adopted values for 




LS 5039 / LS I +61 303 


Jet Lorentz factor (Tint) 


1.02 / 1.25 


Angle between the jet and the observer line of sight (6) 


30° 


Jet velocity (fjot) 


0.2c / 0.6c 


Orbital semi-major axis (a) 


2.6 x 10 12 / 5 x 10 12 cm 


Orbital eccentricity (e) 


0.5 / 0.7 


Distance to the observer (D) 


2.9 / 2 kpc 


Star total luminosity (L s tar) 


1.2 x 10 39 / 2 x 10 38 erg s _1 


Photon flux at the EGRET band (/>ioo Mcv) 


3.3 x 10~ 7 / 8 x 10" 7 photon cm" 2 s" 1 


Photon index at the EGRET band (r) 


~2.2 


Jet radius (R-y) 


10 7 cm 


Injected electron power-law index (p) 


2 / 1.7 


Free parameter 


Adopted values for 




LS 5039 / LS I +61 303 


Magnetic field (B y ) 


1-10 G 


Leptonic jet power (L^ e ) 


10 35 -3 x 10 35 / 10 36 -3 x 10 36 erg s _1 


Maximum electron Lorentz factor (-ySj 1 *) 


10 5 



density variations, implying almost the same effects on the spectrum variability for both cases, 
the EC and the SSC ones. 

These results strengthen the idea that microquasars are 7-ray emitters, being LS 5039 
and LS I +61 303 likely counterparts of the two EGRET sources 3EG J1824-1514 and 
3EG J0241+6103 respectively. However, further observations with the new generation of 7- 
ray instruments (i.e. AGILE, GLAST, HESS, MAGIC, etc..) are mandatory in order to get 
enough angular resolution, flux and temporal sensitivity and energy range width to determine 
properly the association microquasars/EGRET sources and to better constrain the current 
models of emission. 
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Fig. 1 Computed IC spectral photon distribution above 100 MeV for LS 5039 using 
the physical parameter values presented in Table ^ The EGRET data points are 
also shown. The upper limits on undetected EGRET points are plotted with arrows. 
Only stellar radiation density changes due to the orbital eccentricity have been taken 
into account, a) A magnetic field of 1 G has been adopted. Also, the IC spectral 
photon distribution for both the apastron and the periastron passage are shown, b) 
A magnetic field of 10 G has been adopted. 
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Fig. 2 Computed spectral photon distribution above 100 MeV for LS I +61 303 
plotted with the EGRET data points. We have used the physical parameter values 
presented in Table Unlike in Fig.Q] the spectral shape and the variability are the 
same for both B y =l and 10 G. Thus there is only one plot. There are plotted the 
computed I' e , for different orbital radii: a (solid line), the distance at the periastron 
passage (a(l — e), dashed line), and the distance at the apastron passage (a(l + e), 
dotted line) . Now, accretion changes due to orbital eccentricity have been included in 
calculations. 
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